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Abstract 
> 

We use a Stochastic Differential Equation satisfied by Brownian motion taking 

^—i . values in the unit sphere S n —i C W 1 and we obtain a Central Limit Theorem for a 

Q^ | sequence of such Brownian motions. We also generalize the results to the case of 

Cn ■ the n-dimensional Ornstein-Uhlenbeck processes. 
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1 Introduction 

This paper m ay be regarded as an exte nsion to higher dimensions of the 2-dimensional 



CH . study made in I Vakeroudis et al. I (120 111 ). 



CN ■ / (k) \ 

^ ■ We now consider a sequence of Brownian motions 9; ,t>0 , fcel taking values in 

the unit sphere 5 , n _i(c IR n ), all starting from the same point on the sphere. In Section 



[21 we introduce a general representation of 0^ in terms of a Stochastic Differential 
Equation. Using this representation, we describe in detail in Section [3] the limit in law, 
as K — y oo, for the renormalized sum: 



'> 



A' 



- ^—1 

X i ■■ , 



e t 



w 



of these processes, indexed by t > 0, and taking values in W 1 . Of course, one could invoke 
the classical Central Limit Theorem (CLT), at least for the finite dimensional marginals 
of [Zf ,t > 0), as K — y oo. However, with the help of stochastic calculus, there is much 
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more to say about the description of the asymptotics. Finally, in this Section, we remark 
that the CLT can be generalized to Ornstein-Uhlenbeck processes taking values in the 
unit sphere S n -i(c R n ). Three technical points are gathered in an App end ix. 



Further extensions may also be obtained, by following e.g. Ilto I ( 119831 ) or lOchi I ( ]1985[ ) 
and studying for which class of functions /(O) we can obtain a functional CLT such as 
(1T2|) (see below) for /(0*), instead of the unique function /*(0) = @ which we study 
here. 



2 A presentation of Brownian motion in the sphere 



As remarked in IStroock I ( 19711 ) and lYor I ( 19841 ) (eq. (4.j), p. 34), Brownian motion 
(®tjt > 0) in the unit sphere S n -i C M n may be viewed as the solution of a Stochastic 
Differential Equation: 



e t = e + / ^(ej-dB,- 



n — 1 



ds O, . 



In (pQ), Bt = ( B\ ; , i < n J , t > 0, denotes a n-dimensional Brownian motion starting from 



a ^ , while (cr 0,1 (x),x G ^-i) denotes the family of n x n matrices (see e.g. Krylov 
J1980J I). defined by: 

(2) 
(3) 



a°'\x) = (Sij - XiXj) id < n , {x E S„_i) 

and/or characterized by: 

cr 0,1 (x) ■ x = 0, and cr 0,1 (a:) ■ y = y, if y ■ x = 0. 

Note that a 0,1 (x) is symmetric and satisfies: cr°' 1 (a;)cr ' 1 (x) = a 0,1 (x). 
Thus, from ([3]), we deduce that: 

a 0,1 (x)m = m — (m-x)x, mGl", 
(cr 0,1 (a;)m) • (cr ' (x)m') = (m ■ m!) — [m ■ x) {vri ■ x) ,m, m! G M n . 



(4) 
(5) 



3 A Central Limit Theorem for a sequence of Brownian 
motions in the sphere 5 n _i(c R 



Let B^ 1 - 1 , . . . , Q( k \ ... be a sequence of such independent and identically distributed Brow- 
nian motions in the sphere S^-i. We aim for a Central Limit Theorem concerning: 



K 



K 



E ! fc) - E e i 



^) 



(6) 



fe=i 



Adding K equations of the kind of ([T]) term by term, for f ®\ , k < K j, it is immediate 



that: 



Zf = Mf 



n — 1 



ds Zf , 



(7) 



with 



VK k=1 Jo 
Thus, from (0), we obtain: 

Zf = exp 



(8) 



t<_ I ~\Ti-l)t\ f\ m f{n-l)s\ d][/[K 



exp 



o 



(9) 



Now, clearly, the Central Limit Theorem for (Z^), K — > oo, which we are seeking, will 
follow from the limit in law of the martingales (M^,t > 0), as K — > oo. We now state 
both limit results in the following: 

Theorem 3.1 a) The sequence of martingales (M^, t > 0) converges in law, as K — » oo, 
towards: 



M, 



(oo) 



n 



1 - - < 6(0) / VI - e~ ns d/3 s 



77. — 1 



■dS 7 . 



(10) 



where (f3 s , s > 0) is a 1- dimensional BM and (B' s , s > 0) is a (n — 1) — dimensional BM 
taking values in the hyperplane which is orthogonal to 0(0), and £?' zs independent of j3. 
b) Consequently, (Z^,t > 0) converges in law, as K — >• oo, towards: 



7 (oo) 



exp 



n — \ 



i 4 ) I exp V 2 



77, — 1 



dM (oo) . 



'IT 



Proof of Theorem I3.lt Using the Law of Large Numbers, it is not difficult to show 
thatS: 



(M*,t>0) 



(law) 



E 



a'- 



O'H©^)! • d^W = (M™,t> 0) , (12) 



where Q(s) = E a 0,1 {Q s ) is a deterministic matrix, depending on s. The RHS of f)12p 

is a centered Gaussian martingale in W 1 . Before computing the square root involved in 
(TT2"j) . we shall first calculate (see (|2]) for the definition of a 0,1 ): 



E[^\e^)]=(5 M -E[Q l (s)Q,(s)}\ 



j<n " 



(13) 



In order to calculate E 



e«ei J) 



as "naturally" as possible, we consider two generic 



vectors m and m! in M. n , and we compute: 



¥m,m'{t) =E[(m- Q t ) (m' ■ Q t )] . 



(14) 



ft In the Appendix lA.il a more general result, concerning — = L Ylk=i -^» ' dB» is presented, where 
(B^ k \ k = 1, . . . , K) are X independent BMs and {H^ k \ B^) are k < K iid random vectors. 



Using ([[]) and the (special) properties of the matrices {a 0,1 (x)}, we easily deduce from 
Ito's formula, that: 

E [(m ■ t ) {w! ■ e t )] = (m • e ) (m' ■ e ) - (n - 1) / ds E [(m ■ Q s ) (rri ■ 6 S )] 



+ dsE [{a°' l {Q s )m) ■ (a°>\e s )m')] . (15) 

Jo 

Using (J5J, (TT5"j) simplifies as: 

E [(m ■ t ) (m' • 4 )] = (m • 6 ) (m' • O ) - (n - 1) / ds E [(m • 6 S ) (to' • 6 S )] 

+ I ds (m ■ m! — E [(m ■ 6 S ) (m ■ O s )}) 
Jo 

(m ■ 6 ) (m' • 0o ) + (m-m')t-n ds E [(m ■ 6 S ) (w! ■ 6 S )] . (16) 

Consequently, the function ip m>m i{f) = E [(m ■ Oj) (m' • 0t)] is the solution of a first order 
linear differential equation, hence: 



j-t 

n —nt 



E [(m • e t ) (m' • e t )] = e" ni <^ (m • 6 ) (m' • 6 ) + (m • to') / e ns ds ^ . (17) 



Now, taking m = e« and m' = e,-, where (e^; k < n) is the canonical basis of ~R n , the 
matrix Q(s) has elements: 

for i ^ j, (Q(s)) id = -E [Q l (s)Q,(s)] = -0 2 (O)6,(O)e— (18) 

forz = j, (Q( S )).. = 1- £[6,(5)6^)] 

r / p ns _ i 

= l-|(e j (0)) 2 e-" s + e-"N^^— 

- I) + e~ ns ( - - (6,(0)) 2 ) . (19) 



Finally: 



Q( ■ ) = (l--)ld+ e~ ns \-8ii - 6,(0)6,(0) 

n I \ U ' i,j<n 



j<n 



(20) 



= U ~ I (1 - e~ ns )^j Id - e— (6,(0)6,(0)). . 
Using (J20D in the RHS of (02), we obtain: 

M^ = f VW) dB s . (21) 

Jo 



Now y/Q(aj 

A(s) = 



A(s), where (for the explicit calculation, see Appendix IA. 2 j) : 



1 



1 _ _y/l_ e -ns Id+ Wl - - 4/1+ - 

n V n \ \ n — 1 



Vl - e~ ns a ' 1 (9(0)) 



l--< Vl-e-™Id+ W1 + 

U I \ V 71—1 



_ Vl - e-™ ex ' 1 (6(0)) J . (22) 



Thus, (|2Tj) now writes: 



M 



(00) 



A(s) rfS s 







1 - - I / Vl - e- ns dB s + 
« Wo 



1 + 



71—1 



Vl - e" ns 



a^(Q(0))dB s 



(23) 



We remark here that, with /3 S = ©(0) • B s , 

B' s ^B s -e(0)/3 s = a°' 1 (e(0))B s 



(24) 



is a (n — l)-dimensional BM taking values in the hyperplane which is orthogonal to O(0). 

Thus, from (J23J we deduce (fTUj). 

From (jUj), letting AT — >• 00, we obtain ( ITTj) . | 

Moreover, changing the variables s = t — u and using the dominated convergence Theo- 
rem, we have: 



„(oo) (law) I 1 / U- 1 \ /• /72-1 

Z+ = a/1 - — exp — t / exp I — : — -. 1 x 



7? 



x { 9(0)Vl - e-^dPs + \/l + ——^dB' s 



(25) 



s =t-u 1 r ( n - 1 



(law) 



1- 



7? 



exp 



-;v 



p— ra(t— u) 

1 _ e -n(t-«)e(0)d / 9 u + 4 / 1 + — dBj, 

71—1 



t— >oo / 1 1 

1 / exp 



7? 







n 



2 I 



-u] e(0)d/3 u + J 1 - - / exp \-^——u ) (IB 1 .,. 



n 



(26) 



Proposition 3.2 27ie following asymptotic results hold: 
a) 



A 



(00) (law) (30) 



# 



t— >oo 



(27) 



where: 



Z { °° ] = A/1 



n 



n — 1 



w dR, 



b) 



Z\ — exp 



n—\ 



exp 



1 /n-1 \ ,„ i? 

1 exp s dB s — \ 0. 

n V 2 / <->oo 



(2f 



(29) 



Part a) of Proposition 13.21 follows from the previous calculations, using ( QMJ) . In order to 
prove part 6), it suffices to use the expression (125]) and the following Proposition, which 
reinforces the convergence in L 2 result in (129]) . 

Proposition 3.3 As t — > oo, the Gaussian martingales: 

-t 



(Gfe(o),*>o)=e(o)(7 vT 



e- ns e^ s d(3 s - / e^ s dp s ,t> ) , (30) 



and 



(C t) t>0) 



1 + 



-e^W 



n 



eV s dB't> 



(31) 



converge a.s. and in L 2 , and the limit variables are Gaussian, with variances, respectively: 



/ VgX-l+j) n+l I flnrf 



2 a Fi(-i,-U- 



n" 1 n"* 1 — n 



-i 



n— 1 



Proof of Proposition 13.31 

a) The increasing process of the real- valued Gaussian martingale Gj is: 

/ e (n-i)- (yi-e-^-lYds, 

which converges, as t — > oo; thus: 



G 



(0) 



t— >OD 






where the convergence holds both a.s. and in every L p . Of course, the limit variable 
is Gaussian and its variance is given by (we change the variables u = e~ ns and B(a,b) 
denotes the Beta function with arguments a and a**l): 



ds e {n - 1)s VT^ 



n 



du u 2+n (vl 



u 



i) ; 



1 
n 
1 
n 



du u- 2+ n ((1 



-1 + -2 
n 



u 



2B 



2 v / T^+i) 
i 1 3 

-1 + - o 

n 2 



n 



n 



0FT(-l + i) n + l 



/? 



r(i + i) 



n — 1 



*We recall that (T(x),x > 0) denotes the Gamma function and B(a, b) = rta+b) ■ 



To be rigorous, the integral J Q du u~ a {\/l — u — l) , which is well defined for < a < 1, 

can be extended analytically for any complex a with Re (a) < 3. 

b) Likewise, the "increasing process" of the vector-valued Gaussian martingale G' t is: 



y i Hf^-i-^* 



which also converges as t — y oo. The limit variable: 



1 + i—^' - e^ s ) dB' s . 
n — 1 I 

is also Gaussian and, by repeating the previous calculation, we easily compute its variance. 

I 

Proof of Proposition 13.21 : 

From Proposition [3T31 by multiplying both processes by e^ 2 > \ 1 — -, we obtain 



Remark 3.4 (The Ornstein- Uhlenbeck case) 

In fact, for every process satisfying: 

dZ s = dM s + h(\Z s \)Z s ds, (32) 

where (B 4 , t > 0) is a n- dimensional Brownian motion (BM) and h : R+ — >■ R is a bounded 
function, there is a CLT of the kind of Theorem \3.1\ See Appendix \A.S\ for the proof. 

A Appendix 

A.l Generalization for a class of symmetric matrices 

For K independent Brownian motions, and a class of symmetric matrices H ^ such that 
(H^ k \ B^) k < K are iid, we have: 

Mf) = ^-fy H ^.dB^ { ^ fh a dB ai (33) 

with h s a deterministic symmetric positive definite matrix and (B t , t > 0) a n-dimensional 

BM. 

Indeed, using m a generic vector in R™, we have: 

m ■ M? } = m-^/V #(*> ■ dB™ = -L V f (H^m) ■ dB™ 

VKJo k=1 VK k=1 J 

Af(o; f dsE [\H l s ■ m\ 2 ] J , (34) 



(law) 



and for the variance, we have: 

/ dsE [\H] ■ m\ 2 ] =m- I dsE \H]Hl\ m, (35) 

Jo Jo <- 

and 

E [Hi] = h 2 s . (36) 

Remark A.l In our case, we have: H^H^ = Hg. 

A. 2 Square root of Q s 

Q{s) = (l --) Id + e~ ns ( -Sij-Qi^G ^O)) 

V n J \ n J (i,j<n) 

= (l - - J (1 - e- ns ) Id + e~ ns ■ a ' 1 (6(0)) . 

We are searching for a(s) and b(s) such that: 

(a( S )J + 6( S )a ' 1 (6(0))) 2 = Q( S ), 
or equivalently: 

(a(s)) 2 I + 2a(s)b(s)a°> 1 (6(0)) + (b(s)) 2 a ' 1 (6(0)) = Q(s). 
We compare with (137|) and we find: 

(a(s)f =(l-^\(l- e — ) ; 2o( 8 )6(«) + (b(s)) 2 = e~ ns . 

Solving this system of equations, we easily obtain: 



(37) 



(3* 



a(s) = y/1 - JV1 - e~ ns 

,^) = A A T I(v / i-^-v / r^^ 

A. 3 The Ornstein-Uhlenbeck case 

We consider the rz-dimensional 6rnstein-Uhlenbeck (OU) process: 

Z t = z + B t - A f Z s ds, (39) 

Jo 

where (M t ,t > 0) is a 77,-dimensional Brownian motion (BM), Zq G M. n and A > 0. 

Proposition A. 2 The Ornstein-Uhlenbeck (OU) process [Qt,t > 0] in the unit sphere 
S n -i is the solution of the Stochastic Differential Equation 

ef = e z + jTVcef) • d% s - ^1 + A ) J\s ef , (40) 

where (B t ,t > 0) zs a n-dimensional BM. 



Proof of Proposition IA.2b 

We shall study (ft = Sh. We remark that the Jacobi matrix and the Hessian matrix 
associated respectively to the functions $(x) = jfy, (x^O) and g(x) = \x\ are given by: 



*<(*) = 0*°'V) 5 7^ = o* u '>)- 



Hence, using ( 13"9"j) . satisfies the following Stochastic Differential Equation 



<fio + 


Jo Z s 


<^o + 


/ ^s 


00 + 


/ -^s 



n 



1 /•* <fe 



w = vo : / —a^i^-dZs-^j^ J -^-<p 8 111) 

* 0ll fo) ■ rfB s - f (^T§f 2 + ^f^Cft) ) ds 

*' U \ I 5 I I 5 I 

cr ' 1 ^) -Ms- I (p s (^y^^ + A o r0 ' 1 (^) ) ds. (42) 
We can replace the BM B by another BM 1*: 

R* = / 

Jo 

where (Wt,t > 0) is a BM independent from B. 

Thus, f 7 t = f -Mri ■ dM s ,t > 0] and B* are two independent BMs, and from Knight's 



(a°'\<p s ) ■ M s + tr 1,0 (&) • dW s ) , 



theorem (see e.g. iRevuz and Yor I ( 119991 ) and other references therein), 7 is independent 



from the BM (M t , t > 0) obtained by changing the time scale of (J Q yJ-rdB* ) with the 

inverse of f r^n- Finally, (Qf,t > 0) may be obtained from ((pt,t > 0) by making the 
same change of time scale. | 

Corollary A. 3 The angular part of a 2- dimensional Ornstein- Uhlenbeck process is equal 
to the angular part of a planar Brownian motion, considered under the time scale a t = 

e 2Xt -l 
2A - 

Re mark A. 4 For fur ther results concerning the case of a complex-valued OU process, 



see 



Vakeroudis I 1(201 A) 



Proof of Corollary IA.3b 

It follows easily from equation (I3"9"j) for n = 2 by taking the angular part. For the new 
time scale, it suffices to remark that, with < • > denoting the quadratic variation of a 
martingale: 



,"* 2At _ 1 



Remark A. 5 Proposition \A.2\ is easily generalized for every process of the kind: 

dZ s = dM s + h(\Z s \)Z s ds, (43) 

for every bounded function h : M. n — > R. 
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